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At large baryon number density, it is likely that the ground state of QCD is a color-flavor- 
locked phase with a K'^ condensate. The CFL-t-if'^phase is known to support superconducting 
vortex strings, and it has been previously suggested that it may also support vortons, which are 
superconducting vortex rings. We reexamine the question of the stability of vortons, taking into 
account electromagnetic effects, which make leading-order contributions to vorton dynamics but were 
not investigated in previous work. We find that current-carrying and electrically charged vortons 
can be stabilized either by their angular momentum, by Coulomb repulsion, or by a combination of 
both effects. 



I. INTRODUCTION 

The fate of strongly interacting matter at densities 
above nuclear matter remains one of the central open 
problems in the standard model of particle physics. Its 
qualitative behavior is reasonably well-understood, al- 
beit with important inherent uncertainties due to the 
lack of a reliable non-perturbative calculational method 
in QCD. Knowledge accumulated mostly in the last ten 
years indicates the existence of a color superconducting 
phase (for reviews on this vast topic see, for instance 
[1-3]). Color superconductivity arises from the pairing — 
in the BCS sense — of quarks in a certain combination 
of spin, flavor and color; different phases arise depend- 
ing on which quark-quark combinations pair up. At 
asymptotically high densities the pattern of symmetry 
breaking due to the quark pairing in three-flavor QCD 
is known as weak coupling calculations are reliable [4- 
6] in this regime. The symmetry-breaking pairing in- 
volves quarks of all colors and flavors and goes by the 
name of color- flavor- locking (CFL). At lower densities, 
of the kind possibly found in the interiors of astrophys- 
ical compact objects, the problem is more complicated. 
Many proposals for the phases of QCD at these more 
phenomenologically-relevant densities have been made in 
recent years, but a few stand out as appearing more likely 
based on model calculations. Among those few phases, 
the CFL-X*^ phase [7, 8] appears especially likely as its 
existence is suggested by very model independent argu- 
ments at least in some range of densities^. An unusual 
type of solitonic excitation of the CFL-iiT^ phase — the 
"vorton" — is the topic of this paper. 

The low energy degrees of freedom of the CFL phase 
are the Goldstone bosons resulting from the symmetry 
breaking leading to quark pairing. These degrees of free- 
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^ The main loophole in the arguments is the possibility of an 
unexpectedly large strength of (7(l)x-breaking effects at high 
density [9]. 



dom have some similarities to the pseudo-scalar nonet 
found in zero density QCD with one important difference: 
the "kaons" are actually lighter than the "pions" [10-12]. 
The interactions these pseudoscalars exert on one another 
as well as their interaction with electroweak currents are 
largely dictated by symmetry considerations. In fact, the 
low energy effective theory describing them is very similar 
to chiral perturbation theory for QCD at zero density [13]. 
At asymptotically high densities, when the chemical po- 
tential /i is much larger than the quark masses, the quark 
masses can be neglected. As the density is lowered, the 
quark masses, in particular the strange quark mass, be- 
come more important. The effect of the strange quark 
mass can be tracked down to the appearance of a ficti- 
tious strangeness chemical potential for the mesons. In 
particular, the and mesons become lighter 

and, at some critical value of the parameter m'^/^, the 
K° meson becomes massless and condenses, leading to 
the CFL-if° phase [7, 8]. Only small isospin breaking ef- 
fects, like quark masses and electromagnetism, favor the 
condensation of over condensation of K'^ . In the ab- 
sence of a K'-' condensate the promptly condenses. 
Besides the usual vortex solution related to the break- 
ing of the baryon number symmetry U{1)b, the CFL- 
K'^ phase contains another type of vortex related to the 
breaking of the (approximate) symmetry U{l)s^di the 
phase of the field. At the core of this K'^ vortex, 
the K'^ field vanishes and the field may acquire an 
expectation value. 

The condensation of radically changes the prop- 
erties of the vortex. Since is charged, a nonzero 
expectation value for makes the vortex supercon- 
ducting. This mechanism is a specific realization of the 
general mechanism of superconducting strings proposed 
a long time ago in the setting of relativistic field the- 
ory models[14] . Superconducting strings were extensively 
studied in possible realizations in grand-unified theories, 
since they might be formed in the early universe. (For 
some reviews of this vast topic, see for instance [15, 16]). 
Kaplan and Reddy pointed out that the CFL-i^° phase 
contains all of the ingredients necessary for the formation 
of superconducting strings [17]. More importantly, there 
is a natural mechanism for the stabilization of supercon- 
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ducting vortex string loops, which are often referred to 
as "vortons" . If there is a supercurrent flowing through 
the loop, there is an associated energy proportional to 
~ RI^, where / is the current and R the loop radius. As 
the radius of the loop changes, the current is not con- 
stant; it is actually proportional to 1/R. Naively, this is 
because supercurrents are related to the winding of the 
phase of around the loop, a topological invariant. As 
we will see here, the argument becomes more subtle once 
gauge invariance is taken into account, but the behavior 
/ ^ R~^ is maintained by the conservation of angular 
momentum. Against the contribution of a supercurrent 
to the energy of the loop, proportional to PR ~ l/R, 
one must balance the energy of the vortex due to its ten- 
sion, proportional to R. One then expects there to be 
some finite value of R that minimizes the energy, and is 
the stable radius of the vorton. Buckley et al. pointed 
out that this mechanism did not guarantee the stability 
of quark matter vortons[18]. The reason is that the ra- 
dius of the vorton should be larger that the thickness of 
the string. One may think that one can always arrange 
for the radius to be larger by making a larger current to 
flow through vorton. However, a large current has the 
effect of quenching the superconducting condensate and 
destroying the vorton, a phenomenon well known in the 
context of cosmic superconducting strings[15]. On the 
other hand, the presence of electric charge in the loop has 
an anti-qucnching effect that also resists the shrinking of 
R and may stabilize the vorton. Only a more detailed 
study, including numerical estimates, can tell whether 
vortons are likely to exist in the CFL-X° phase. 

There are some important differences between vortons 
in the CFL-iiT^ phase and superconducting cosmic strings 
besides the trivial difference in energy scales. First, the 
CFL-K'^ phase shields electric fields very weakly, at least 
at the small temperatures relevant to neutron stars, in 
contrast to the strong shielding at cosmologically high 
temperatures usually assumed in the study of cosmic 
strings. Second, the low energy dynamics of the CFL and 
CFL-if" phases are described by a non-linear cr-model. 
As usual, the sigma model is written in terms of a com- 
pact variable S ~ exp(i7r//7r), so that there is maximum 
value ^ /tt for the meson fields tt. Since the meson fields 
acquire expectation values of roughly the same magni- 
tude as this maximum value, it is unclear whether a lin- 
ear approximation to the non-linear cr-model, which was 
used in previous studies of vortons in the CFL-K^ phase, 
will suffice for the stability analysis. It is the purpose of 
this paper to present the analysis in the full non-linear 
sigma model. Additionally, in previous studies, electro- 
magnetic effects were not fully explored. Here, we will 
see that these effects make leading-order contributions 
to the stability of the vorton. Of course, this should not 
be entirely surprising: the effects of isospin breaking due 
the quark mass splittings, which are so crucial for the 
very existence of superconducting strings in the CFL-X° 
phase, are known to be roughly of the same scale as the 
effects of electromagnetic interactions. 



II. CF'L+K° PHASE 

When quark masses can be neglected as compared to 
the chemical potential /i, the high-density ground state 
of 3-fiavor QCD is believed to be the color-flavor locked 
(CFL) phase. This phase is characterized^ by a conden- 
sate of the form[19] 

(qhCql,) - - (qhCQR,) = m'A e'^'^e^.z (1) 

where q are the quark fields which carry color 
a,b = 1,2,3, fiavor i,j — 1,2,3, and helicity R (right- 
handed), L (left-handed) indices. The size of the con- 
densate A is referred to as the 'gap'. Such a condensate 
breaks the symmetry from SU{3)c x SU{3) l x SU{i)R x 
U{1)b X U{1)a SU{3)c+L+R X Z2 X Z2. The first 
Z2 factor is the remnant of U{1)b that leaves the CFL 
condensate invariant, while the second Z2 factor is the 
remnant of the approximate U{1)a symmetry of QCD 
at high density. There are 18 = 2*8 + 2 broken gener- 
ators, but 8 of the would-be Nambu-Goldstone mesons 
are eaten by the color gauge fields. As a result, the glu- 
ons get masses of order g^^ (g^ is the strong coupling 
constant) [10]. The CFL condensate breaks the vacuum 
U{1)em C SU{S)l+fi. gauge symmetry, but a linear com- 
bination of Ag color generator with the U{1)em genera- 
tor annihilates the condensate, generating an unbroken 
U{1)q gauge symmetry. We will henceforth refer to the 
resulting U{1)q dynamics as electromagnetism for sim- 
plicity'^ . 

The remaining ten broken generators are associated 
with physical Nambu-Goldstone modes. Two of these are 
associated with the superfluidity resulting from U{1)b 
and U{1)a breaking, and will not be important in what 
follows. The other eight form a multiplet transforming 
under SU{3)c+l+r the same way the zero density pseu- 
doscalar octet transform under flavor SU{3). At energies 
small compared to A, the low-energy dynamics of the 
CFL phase is dominated by these physical NG modes. 
The dynamics of the NG "mesons" that arise in the CFL 
phase are described by an effective field theory [11-13] 
which bears close resemblance to SU{3)f chiral pertur- 
bation theory. 

The effective theory can be written in terms of a chiral 
field S, which transforms as E — g^Yign under SU (3) l x 
SU{3)ii- The eight CFL mesons tt", a = 1, . . . , 8 appear 
as fluctuations of the E field around its vacuum value Eq, 



^ The CFL condensate shown above is not gauge invariant. There 
are gauge-invariant order parameters characterizing the CFL 
phase [3], but much as in the electroweak sector of the Standard 
Model, for many purposes it is convenient to discuss the physics 
using non-gauge-invariant language. 

^ In fact, at high density, it turns out that the Ag part of the 
generator of it U{1)q is small when oem is small[20]. Also, the 
differences between oq and aeM, Eind the dielectric constant 
eCFL in the CFL phase and evacuum are both suppressed by the 
smallness of oiem ~ 1/137. 
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and are packaged inside E as E = ^ ^■'^Eq. Here A" 
are the standard Gell-Mann matrices with Tr (A'A^) = 
26^^ , and / is the equivalent of the pion decay constant. 

The lowest-order effective Lagrangian describing the 
CFL+K" phase is [7, 10-13] 

- ^Tr (VoEt . VoE - v^D,^^ ■ D,E) 

+ 2ADet (M) Tr (M^^E + h.c.) - ^Ff.^F'"', (2) 



where 



VoE = L»oE + i 



2pf 



(3) 



Pf is the Fermi momentum (which we take to 
be /i), Af^ is the electromagnetic potential, and 
Q = |diag (2, — 1, — 1) is the quark charge matrix un- 
der the U{1)q symmetry. F^^ is the electromag- 
netic field strength tensor, while M is the mass matrix 
M = diag (m„, m^, m^) « diag (2, 5, 95) MeV. 

The dependence of the mass term is different from 
the mass term encountered in zero-density chiral per- 
turbation theory, which is linear in M. The reason for 
this can be traced to the fact that at high densities, 
the dominant source of chiral symmetry breaking is the 
CFL condensate rather than the usual chiral conden- 
sate. The CFL condensate preserves left-handed and 
right-handed quark number modulo two, which forbids 
mass terms with odd powers of M. At lower densities, 
one expects that a {(Ib^l) condensate will become non- 
negligible, leading to the appearance of 0{M) terms in 
the effective theory. Here we assume that the density is 
large enough that the 0{M) terms are negligible. The 
0{M^) term in eq. (2) is not the most general ©(M^) 
term consistent with the symmetries. As shown in [7], 
however, the other two possible terms are suppressed by 
a power of A//x compared the one in eq. (2). 

In contrast with chiral perturbation theory, the low- 
energy constants appearing in the CFL effective theory 
are calculable in the /i — > oo limit. At very large /x, one 
finds [10] [12] 



21 - 8 In 2 ^ 
18 2^' 



.^ = 1, and A= AA^ (4) 



where /i is the chemical potential and A is the gap. 
The equation of motion for E is 



P 



(VoVqE - v'^D.DiY) = 2AT)et (M) 



-It 



(5) 



The equations of motion for and A* are 
^ -^Tr (V"E^ • [Q, E] - h.c. 



Tr m^E 



h.c. 



(6) 



It is straightforward to read off expressions for the elec- 
trical charge and current. We note that there is a sub- 
tlety: a field that has no time dependence still carries 
charge. In fact, the time-dependence of the fields is not a 
gauge- invariant concept. Only the combination of time- 
variation, charge chemical potential and gauge potential 
is gauge-invariant. 

For values of M small compared to /z, the vacuum value 
of E is simply the unit matrix and the meson mass spec- 
trum can be read off from an expansion around Eq = 1. 
This is the CFL phase. In this phase, the charged pions, 
charged kaons and neutral kaons have masses [7] 



mK± 




(7) 




which arise from the O (M^) mass term in the effec- 
tive Lagrangian. The neutral pion and eta mesons have 
masses of a similar form, but as they lie on the diagonal 
when other fields vanish they commute with the MM''' 
matrix in the covariant timelike derivative and thus are 
never enticed to condense. 

As one considers larger values of quark masses, the 
\MM^ /2pF, E] terms in the Vq covariant derivatives be- 
come non-negligible, and act as a chemical potential 
for the mesons, effectively shifting the kaon masses by 
±.m1/2pF'^- For values of rUg larger than a critical value 

of order ~ rn^J^lS?/^, the neutral kaons mass-squared 
would be negative. Instead, the neutral kaons condense 
and the field E picks an expectation value Eg 7^ 1. Since 
K'^ is a neutral boson, this is a superfluid phase. 

The is also driven to condense by the strangeness 
chemical potential (assuming the the opposite charges re- 
quired by charge neutrality are energetically cheap) . The 
question then arises as which field, or actually 
condenses. Two factors favor the condensation of . 
First, the small quark mass difference — niu makes 
the neutral kaons slightly lighter than the charged ones. 
Second, a K'^ condensate carries electric charge and, in 
order to maintain charge neutrality, electrons have to 
be present. The electrons have an energetic cost due 
to their mass and kinetic energy up to the Fermi level. 
For this reason the condensation of is favored. There 
is, however, a counteracting effect: fluctuations of the 
electromagnetic field lower the energy of charged con- 
densates. A full calculation of this effect in the CFL-if" 
phase requires the knowledge of the coefficient of an ad- 
ditional term in the effective Lagrangian for the mesons 



The shift in the pion masses due to the "strangeness chemical 
potential" term is small, of order ±(m'i — m^)/2pp[7]. 



4 



contributing to the electromagnetic correction to the K'^ 
mass [21]. For most reasonable values of this coefRcient 
it is the field that condenses and we will assume in 
this paper that that is the case. 

With all other fields vanishing, the vacuum expectation 
value of condensate is given by [7] 

where mj^o is the neutral kaon mass coming from the 
0(WP) mass term, while iiij = (mf — m^)/(2yu) is the 
effective chemical potential coming from the presence of 
M matrices in the time-like covariant derivatives. This 
expression determines \K^\ if Msd > ^k": if l^sd < n^K'^ 
there is no if ° condensate. 



III. VORTONS 

The CFL+if" superfluid supports global vortex strings 
as a consequence of the breaking of if" number. At the 
center of such vortices the vacuum expectation value of 
if" vanishes due to the single-valued nature of the field. 
As we will see shortly, in the CFL-fif" phase, the cores 
of vortices are large enough to be describable within the 
mesonic EFT. The vanishing of if" in the cores of vor- 
tex strings provides room for the next lightest boson, the 
if + meson, to condense. Thus, the cores of the if" vor- 
tex strings may contain a charged condensate, which is a 
superconductor. These vortices are therefore supercon- 
ducting vortex strings [14], and may carry electrical cur- 
rent and charge. Since the pions are significantly heavier 
that the kaons we will work in the approximation where 
the pion fields are turned off. We will keep the U{1)q 
gauge fields and the if",if+ fields, and include all of 
their non-linear interactions. Previous studies of vortex 
strings and vortons in the CFL-if " phase did not take 
into account electromagnetic effects on the dynamics of 
the strings, and worked only to quartic order in the me- 
son fields for simplicity. 

To describe a single long if" vortex in the z direction 
the if" field takes the form 

if" = if"(r)e^'*. (9) 

where (j) is the axial coordinate, r the axial distance, and 
I is the vorticity. For definiteness, we henceforth take 
I — 1. Since the vortex may have an overall electric 
charge per unit length Q± and can carry an electrical 
current i, an appropriate general anzatz for the if + field 
takes the form 

if+ =if+(r)e*("*-'=^). (10) 

where k is the wavenumber along the z direction, t is time 
and the frequency w contributes to the charge density. 
Since a if"*" field of the form above carries charge and 
current, it generates an electromagnetic field. The energy 




FIG. 1. A vorton with arrows indicating the direction of phase 
change of the if" (0) and (kz) fields. 



of a vortex configuration is higher that the ground state 
so there is an energetic cost in creating a vortex. The 
excess of the vortex energy per unit of length is the string 
tension. The string tension is actually infrared divergent, 
a well known feature of global vortices. A finite energy 
is obtained if the vortex exist in a finite container or if 
the vortex forms a closed loop. We are interested in the 
second case, and we will think of the z direction as the 
direction along the loop. 

We now consider vortons, which are closed circular vor- 
tex loops. Generically, the string tension of a vortex loop 
tends to make such closed strings smaller. On the other 
hand, if the vorton carries current, the kinetic energy of 
the if"*" field due to its (covariant) derivative along the 
z-direction grows as the vorton shrinks. In addition, the 
presence of charge will help stabilize the vorton due to 
Coulomb repulsion and, indirectly, by helping stabilize 
the if'^ condensate supporting the current. There are 
also magnetic effects to consider as the current on op- 
posite sides of the loop repel each other. Finally, there 
is a contribution to the energy of a vorton due to weak 
interactions that, for a very large vorton, can be appre- 
ciable. The C/(is(l) symmetry related to if" number, as 
well as the Uds{i) symmetry related to if + number, is 
explicitly broken by weak interactions. The consequence 
for the structure of a if" vortex is that it is energeti- 
cally favorable for the if" phase to vary with the angle 
around the vortex in a inhomogeneous way. In fact, it 
is favorable for the phase to change little with the angle 
except in a narrow range of angles where the phase shifts 
somewhat abruptly by 27r [22]. The consequence is the 
formation of a domain wall attached to the vortex carry- 
ing some amount of energy. For a closed vortex forming 
a vorton this wall is a membrane stretching across the 
vorton (which justifies the name "drum vortons" some- 
times used to describe these objects). The tension of the 
domain wall is, however, a weak interaction effect and 
consequently very small. It is only important for very 
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large vortons. As we will see, for the vortons considered 
in this paper the membrane is a negligible effect. 

All the contributions discussed above are important in 
different regions of parameter space and should be taken 
into account. Obviously, finding numerical solutions of 
the full set of coupled non-linear partial differential equa- 
tions describing a closed loop is quite challenging. They 
are hard to find even in the case of a straight vortex. 
Fortunately, in order to get a semi-quantitative under- 
standing of the physics, one does not have to solve the 
full 3D problem. 

Before tackling the problem of determining the me- 
son and electromagnetic field solutions, let us discuss the 
conserved quantities constraining the time evolution of 
vortons. The electric charge of a vorton can only change 
through the absorption of an electron and the conversion 
of a into a . The rate of this weak process is 
suppressed by the fact that, as we will find out, typical 
vortons occupy a volume much smaller than the size of 
the electron clouds around them. There will be a time 
scale, however, where this process will be equilibrated. 
For times larger than that scale the electric charge of 
the vorton will be driven to zero, the condensate 
will also vanish and the vorton will collapse. Our cur- 
rent study is valid for shorter time scales. We leave as a 
future project to estimate the time scale for charge equi- 
libration. A similar statement can be made about the 
angular momentum of the vorton. There are processes 
(quasi-particle and/or photon emission) that can change 
the angular momentum of the vorton. We will consider 
a time scale shorter than that, where the angular mo- 
mentum of the vorton can be considered to be constant. 
There is one more fixed quantity: the winding number of 
the K'^ field along the vorton. 

If the field is to be periodic in z while remaining 
single-valued, k in eq. (10) must be quantized so that 

211 N = I dzk (11) 
Jo 

where N is an integer and R is the vorton radius, so 
that k — N/R. The winding number N cannot change 
smoothly, since the topology forces it to be an integer; it 
can only change by tunneling. One expects the tunneling 
rate to be very small, suppressed by the exponential of 
the energy it costs to set the condensate to zero ev- 
erywhere along the vorton, which is what is necessary to 
change N. Therefore, one expects iV to be a conserved 
quantity. This intuition is correct. However, as we now 
discuss, specifying N does not fix the properties of the 
vorton. 

The quantization of k was already discussed in 
Ref. [18], which did not treat electromagnetic effects. 
These effects are somewhat subtle. The time-dependent 
phase of the field in eq. (10) can be changed by U{1)q 
gauge transformations, and is thus is not an observable. 
However N = kR remains a physical observable, since it 
cannot be changed by periodic gauge transformations. In 



fact, given a loop of superconducting wire, the quantiza- 
tion of the fluxoid connects N to the magnetic flux $b 
that threads the superconducting loop and the current 
density at the loop's center by [23] 

27rjV ^e^B- idz ,, J ^ . (12) 

/ /^Tr(|[Q,S]P) 

In most familiar condensed matter systems, the super- 
conducting sample is large compared to the London pen- 
etration length 1/m^ so the current that flows through 
the center of the loop is exponentially small. Then the 
integral term in eq. (12) is negligible and it then appears 
as though flux itself is quantized. While qualitatively the 
flux and fluxoid are different entities, quantitatively they 
are very similar in this 'thick-wire' case. In our case, we 
will find that the penetration length is comparable to the 
thickness of the wire; the vortex is in neither the "thick" 
or "thin" wire regime.^ Thus, for us the distinction be- 
tween the fiuxoid and the flux is quite significant. 

Next, even though in principle and N/R (or more 
properly $s and N) are independent physical observ- 
ables, they do not enter calculations separately. Instead 
they appear in just such a combination so that they al- 
ways give the current density j^; this can be traced to 
the fact that the contributions of and N/R both orig- 
inate in the U{1)q covariant derivative — dz + A^. 
This poses a puzzle, because although N is conserved, it 
does not determine the current, and we cannot deduce 
how the current should change as the radius of the vor- 
ton changes. If we increase A'^ by 1 and the magnetic flux 
by 27r/e, the resulting vorton will have the same current, 
and so the Lagrangian will not differ for these two differ- 
ent vortons. 

In the ungauged case the dependence of the current 
density j on the vorton radius R is quite straightforward: 
j k ^ N/R. This scaling was, in fact, used to argue for 
the stability of the ungauged vorton. Once we take the 
gauge fields into account, it is not immediately obvious 
that this stabilizing effect will persist. Conceivably, the 
gauge field might adjust in such a way as to cancel the 
R~^ behavior. As will be made explicit in section IV, 
a stabilizing effect of this sort survives, and shows up 
through the requirement of angular momentum conser- 
vation. 

We have argued that in order to describe a vorton, we 
should specify both N and the magnetic flux — these 
quantities together give the current. If we specify both 
the current and the charge, we can compute the angu- 
lar momentum. Instead, we will just specify the angular 
momentum and the charge. Thus, when we specify the 
angular momentum, we are actually discussing a whole 
class of vortons whose magnetic flux and winding num- 
ber conspire to give the specified value of the angular 



Cosmic strings can typically be analyzed in the much simpler 
"thin" wire approximation [14] 
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FIG. 2. A cross section of the three regions of integration 
described in the text. The inner torus (orange online) is the 
region where the energy is approximated by 2'kR times the 
tension. Outside the sphere (green online) the fields are ap- 
proximated by the multipole expansion. The apple-core shape 
(blue online) is neglected. 



momentum. While vortons supporting different amounts 
of magnetic flux and winding number are in principle 
physically distinct objects, the stability and equilibrium 
properties of vortons related by having the same angular 
momentum must be the same. 

Thus, we are interested in minimizing the action at a 
fixed charge and angular momentum. One way to ac- 
complish this is to use Lagrange multipliers v and fl to 
define 

F = i^Q + QJ - Leff. (13) 

Extremizing F is equivalent to solving the equations of 
motion arising from the Lagrangian L while guaranteeing 
that the field configuration supports the correct amount 
of charge and angular momentum. We will need to elimi- 
nate these Lagrange multipliers appropriately in favor of 
given values of the charge Q and angular momentum J. 

The calculation of F by integrating its density over 
space is very complicated due to the geometry of the 
vorton. We will obtain an approximate solution to this 
problem by dividing the space into three regions. One 
is given by a torus of radius R around the vorton (see 
Figure 2). The second is the region outside a sphere 
of radius 2R centered on the vorton center. The third 
region is the small "apple-core" shaped region between 
the torus and the sphere. By picking the sphere's radius 
to be 2R, we ensure that it just touches the torus but 
that they do not overlap. The integral over the torus 



will be approximated by 27ri? times the quantity 
Jo 

fR 

= vQl_ / (fTl_ Ceff (14) 

Jo 

where Q and are the electric charge and angular mo- 
mentum densities, and Q± and J'± are the corresponding 
linear densities. The quantity T± plays the role of the 
string tension when charge and angular momentum are 
fixed. We will compute from the definition above by 
taking the fields to be the one in a straight vortex. This 
approximation becomes exact as the vorton loop radius 
R is much larger than the vortex thickness S. The in- 
tegration over the region outside the sphere is simple. 
The charge density will vanish in the ansatz for the kaon 
fields we will adopt (and it is very small in the exact solu- 
tion). The only contributions to the angular momentum 
will come from long-distance electric and magnetic fields, 
treated in the multipole expansion. The "apple-core"- 
shaped region between the torus and the sphere will be 
neglected. In order to proceed we will now consider the 
calculation of J-± for a straight vortex. 

III.l. Straight superconducting vortices 

In the presence of a single, infinitely long straight ff" 
vortex in the z direction, the K''^ field takes the form 

K" = K^{r)e'^'^. (15) 

where (j) is the axial coordinate, and I is the vorticity. 
For definiteness, we henceforth take / — 1. Since the 
vortex may have an overall electric charge per unit length 
Q_L and can carry an electrical current /, an appropriate 
general anzatz for the field takes the form 

K+ = K+ (r)e'^'^'-^''\ (16) 

where k is the wavenumber along the z direction, t 
is time and the frequency u contributes to the charge 
density. Since a K'^ field of the form above car- 
ries charge and current, it generates an electromagnetic 
field. The equations of motions following from eq. (2) 
then give four coupled non-linear ODEs that determine 
K\r),K+{r),At{r),A,{r). 

Unfortunately, these equations cannot be solved an- 
alytically. Nevertheless, we can get some insight into 
the behavior of the solutions by posing a simple form 
for the radial profiles of the kaon fields and calculating 
the electromagnetic response to these background pro- 
files. On general grounds, one expects that the mag- 
nitude of the condensate in a single-vortex solution 
will interpolate between the limiting values K^{Q) — 
and K'^ir — >■ oo) ~ |iir°| which takes the value given by 
eq. (8). At the same time, the K'^ condensate profile will 
go from some fixed value = \K^\ at r = to zero 
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at large r. The simplest profiles for the kaon fields that 
incorporate these general features is 

K°{r) = \K'^\e{r-5), K+ {r) = \K+\e{5 - r), (17) 

where we will refer to 5 as the string thickness and treat 
it as a variational parameter alongside the parameter 
|iir^|. Obviously, these simple step- function ansatze will 
not solve the kaon equations of motion, but we expect 
them to capture the qualitative features of the full solu- 
tions. While somewhat simplistic, these square profiles 
provide the ability to solve the electromagnetic equations 
of motion and understand their back reaction on the kaon 
fields analytically. 

Given the above ansatze we can solve the electromag- 
netic equations of motion in closed form. In keeping with 
our anzatz, we do not allow the electromagnetic fields to 
backreact on the shapes of the kaon fields, but we do 
allow the electromagnetic effects to change the size of 
the charged condensate at the center of the vortex. In 
solving for the profiles of the gauge fields, it is useful 
to note that k,uj and At,Az and ^su are not separately 
gauge- invariant. They appear in the equations of motion 
through the gauge-invariant combinations 

Co{r) = uj ~ eAa{r) - ^lsu k{r) ^ k + eA^ir). (18) 

Alone, the time dependence lu, the size of the gauge field 
and the quantity fisu are meaningless: explicit time 
dependence and the chemical potential may be absorbed 
into an overall shift in Aq and only grouped into ui do 
they take on a gauge-invariant meaning. An analogous 
statement may be made about k and Az- In the static 
solution we are looking for we may assume that the gauge 
fields are time independent. Then, in Coulomb gauge and 
with the ansatze in eqs. (15), (16) and (17), the gauge 
field equations of motion (6) reduce to 

V^cj = Cbm^ e{S ~ r) 

V^k = kv^m'^e{S -r), (19) 

where = f^e^ sm^{\K^ \ V^/ f) is the mass of the pho- 
ton inside the vortex core (outside of the vortex core the 
photon is massless). We also see that the photon mass is 
bounded by — e^/^, which is an effect that would be 
missed if one expanded E in the kaon fields. The charge 
and current densities may be read off directly from the 
right side of these equations. Notice that the charge and 
current depend on the value of the electromagnetic fields 
themselves. 

The electromagnetic field equations of motion, by 
gauge invariance, are homogeneous equations of u and 
k respectively. In contrast, when electromagnetism is 
sourced by spin-half fermions, the equations of motion 
are not homogeneous — the electric charge is proportional 
to the fermion number alone. Instead, we find that a 
nonzero field only carries charge if the gauge field is 
nonzero. This implies that if the charge density vanishes 



at the wire's core there cannot be electric charge any- 
where. This is apparent in the explicit solutions to eq. 
(19). An analogous statement holds for current density. 
The solutions to the equations of motion are 

uj{r) = ujQ It){m^r) r < 5 

= ujQ (^lQ{m~^S) + ■m^SIi{m^S) log (^-^^ r>S 

(20) 

k{r) = ko loivnijr) r < 5 

^ ko (^Io{vmyS) + vm^SIi{vmyd)log (^-^^ r>5 

(21) 

where Iq and Ii are the O*'' and 1"* modified Bessel func- 
tion and Wo = uj{r = 0) and k^ = k(r = 0). The values of 
(Iiq and fcg are fixed by the boundary conditions. We can 
relate them to the charge and current flowing through the 
string. In fact, using the known large-r behavior of the 
electromagnetic fields due to an infinitely-long charged 
current-carrying wire, we find that the charge per unit 
length Q_L and the total current / carried by the vortex 
are given by 

Qj^ = H LOQm-^SIiirn^d) — LUQmjSfi{m-yS), (22) 

2tt ~ 

I = kQvm~f5Ii{vm^5) = —kovm-yS fi{vmyS) 

(23) 

We emphasize that, due to gauge invariance, the solu- 
tions only depend on w, fisu , k , Aq and A^ through the 
combinations co and k. A change in, say, w is absorbed 
into a change in Aq in order to keep uj the same. The ap- 
parent difference in sign between the charge and current 
is due to the difference in sign of the gauge fields when 
defining w and k. 

We now calculate the Lagrangian per unit of vortex 
length C± for the gauge fields in eqs. (20,21). The inte- 
gral in the radial direction is infrared divergent due to the 
logarithmic behavior of the gauge field at large distances. 
As discussed above, we will only need to integrate the La- 
grangian up to a radial distance R. The Lagrangian per 
unit length £± is then a function of {l^sT"''! ,ojo,ko,5,R} 
defined by 

C± = / d^rj_ £eff. (24) 
Jo 

A calculation with the ansatz in eq. (17) and the solutions 
in eqs. (20, 21) gives a somewhat unwieldy but analytic 
expression 

C± — CjQim^S fj.^m^S) — k^vm^S f^{vm^5) + Afr^, (25) 
where 

hix) = ^ {lo{x)h{x)+xI^{x)logiR/S)) , (26) 
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and 



+ 47r(i?- - (5 ) mams + mu{md + TOs) 




2A 



{l\^\og{Rl6)) 




(27) 



The dependence on is hidden inside the photon mass 
™7 = e/ sin(-\/2if "*"//)■ The first two contributions to 
/5 come from the integral of the mass term in the La- 
grangian. The second two come from the kinetic 
terms. Incidentally, we mention that there are electro- 
magnetic one-loop radiative corrections contributing to 
the energy split between a and a K'^ condensate[21]. 
The magnitude of this effect is comparable to the one 
given by the difference in quark masses included above. 
Unfortunately, the calculation of its precise value de- 
pends on the value of a counterterm in C^f / not yet com- 
puted even in QCD perturbation theory. For reasonable 
values of this counterterm this contribution is not suffi- 
cient to revert the roles played by K'^ and In what 
follows, we will disregard this contribution. 

Extremizing £_l with wq, fixed leads to an equation 
of motion for \K^\. Depending on the values of wcfco 
two kinds of solutions are possible. One class of solu- 
tions has a nonzero vacuum expectation value and is 
associated with the existence of a charged condensate at 
the core of the vortex which breaks the U{1)q symme- 
try. The other class of solutions has a quenched conden- 
sate, satisfying sin(-\/2 |A'+| //) = 0, and does not break 
U{\)q. These other insulating solutions correspond to 
E = diag (—1,1,-1), which is diagonal and so commutes 
with Q. This implies that the photon mass vanishes, as 
can be seen from eq. (6). The existence of this pecu- 
liar quenched solution is a direct result of the kaon mag- 
nitude being a compact variable, and is not seen upon 
linearization. Despite this odd property of the full non- 
linear theory, it remains the case that a large wavenum- 
ber generically causes quenching. This is because one can 
show that if the preferred value of l-^f"*"! //is near 0, 
increasing fco tends to push the condensate smaller; if 
the preferred value is near tt, increasing fcp tends to push 
y/2\K+\/f toward tt. 

Part of the dependence of the solutions on / and 
Q_L is easy to understand and was already anticipated 
in previous work which used a linearized version of 
eq. (2) [18, 24]. A large |fco|, which enables a large 
current, simultaneously decreases the magnitude of the 
charged condensate and eventually leads to quenching, 
while a large which helps increase Qj^, helps stabi- 
lizes the condensate. 



If our main interest were the straight vortex we could 
now minimize >C_l in relation to \K^\ and 5 and trade the 
dependence on loq and by Qj^ and /. That would tell 
us whether or not the condensate is quenched inside the 
vortex and how thick the vortex is in a situation where 
the charge density and current are held fixed. This is, 
however, a very different problem from the one we aim 
to solve. In a vorton with changing radius it is the total 
charge and angular momentum that is held fixed, not 
the charge density and current. For this reason, we will 
proceed to calculate the remaining terms contributing to 

The charge per unit length is given by the expression 
in eq. (22). The angular momentum per unit of length 
requires a little more work. The total momentum density 
for the field configuration is 

f 



j.Ot ^ J rj.^. (vD'J:^ ■ V"E + h.c.) + ExB 



(28) 



We are interested in the momentum that fiows along the 
vortex. With our ansatze, we find 



T'^=C.,~k,vf sin' nJ^^ 



(29) 



{Io{mjr)Io{m^vr) + Ii{m~fr)Ii(m-yVr)) r < 6 
S' 

-^Ii{m^S)Ii{m^vS) r > S, 

The contribution from outside of the vortex is provided 
entirely by the electromagnetic Poynting vector. The to- 
tal angular momentum per unit length is given by 

Jj_ = y dV_L ^Jr' + rl + 2Rr^ cos </> (30) 

R 

dV_L R 



Rmy6koU]of2, (31) 



with 



2Trv 



e2(l + v) 



{lQ{vm^S)Ii{m^S) + Ii){m^5)Ii{vm^5) 



+ (1 + v)m^5Ii[m^5)Ii{vm^5) log{R/5)), (32) 

where we have inserted the lever-arm from the vorton's 
center to the vortex's center, R, to approximate the an- 
gular momentum per unit length. Corrections due to 
different portions of the vortex being different distances 
from the vorton's center are negligible in the limit R^ S 
since most of the angular momentum is carried inside of 
the vorton (r < S) and not by the electromagnetic field 
outside of the vorton as we numerically verified. 



III. 2. Far fields 

Outside a sphere of radius 2R the main contribution 
to L and J come from the electric and magnetic fields. 
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We will approximate the fields in this region by their 
multipole expansion: the electric field by the monopole 
contribution and the magnetic field by its dipole form. 
The fields are 



E. 



Bout = (3(m • f)f - rn) 



(33) 
(34) 



Q — 2ttRQ± is the total vorton charge, and m = ttR^ I z 
is the magnetic moment of the vorton, assuming it lies in 
the xy— plane. 

Their contribution to the total Lagrangian outside of 
the sphere is 



in 



outside 
oo 



2R 



Q2 ^2 



167ri? 967ri?3 
= {2ttR 6m.,) ((1.2 /7 - v~kl fs) (35) 

where we have grouped various factors to appear similar 
to (25) when we evaluate 2TrR£±, so that 



TT^ m^(5 2 



2 e2 
48 e2 



/i(wm^(5)2 



(36) 



The difference between /y and /g, aside from the factors 
of V that are associated with ko, is the dimensionless 
factor that arises from the fact that the magnetic field is 
a dipole and not a monopole like the electric field. 

Further contributions from multipole terms will be 
down by at least factors of 4, which correspond to the 



factor of 2 in front of the cutoff 



= 2i? as well as 



geometrical factors. The higher electric and magnetic 
multipole moment contributions will have the effect of 
changing the coefficients of electric monopole and mag- 
netic dipole terms at C (1), but will not change the qual- 
itative dependence on R, Q, and J shown above. It may 
seem strange that the contributions of the electric and 
magnetic energy have opposite signs as both should be 
repulsive and help vorton stabilization, and thus presum- 
ably should have the same sign when they enter F. This 
impression is an artifact of the way we set up our calcu- 
lation up to now. After the Lagrange multipliers v and 
il are eliminated in favor of the total charge both electric 
and magnetic energy will have the same repulsive effect. 

The angular momentum also has a contribution coming 
from the far fields as the Poynting vector S — Eout x Bout 
field lines forms closed circles around the vorton. This 



contribution can be calculated as 
Jout = d^r fx S 

J outside 



r^dr dfl fx | - ( 
2R V (47r)2r5 

Qm 



with 



12nR 

(2TrR)RmjSkoUJo fe, 



(37) 



(38) 



arranged to be easily added to the contribution to the 
angular momentum from 2ttRJ^±. 



IV. VORTON STABILITY 

We can now collect the results in eqs. (22,25,30,35) and 
(37) and obtain the function to be minimized F by adding 
the contributions coming from the three regions of space 
(see Figure 2). F is a function of {l-ftT'*'! ,u}o,ko,6, R, fl} 
but not Q or J — these quantities are expressed wholly 
terms of those variables as well. We find 

F — V (27ri?a;om^(5/i(m-y5)) -f n{2TrR'^mj6koojo){f2 + /e) 

- 27ri?(5m.^(I)o (/3(m^(5) + fj) 
+ 2nR5myvkl {fsivm^^S) + /§) 

- 2nRAh (39) 

We now minimize F in relation to ujq and fcg. This gives 
us two relations that we can use to eliminate v and f2: 



<^oR f2 + h 



Plugging these values of ly and il back into F we find 



F = 2TTRm.,S 



iMm-,S) + h) + P ifiivm^S) + fs)) 

(41) 



2nR A/5 



Having eliminated the Lagrange multipliers from F, we 
can trade ujq and fco for the conserved quantities. The 
quantities we are given by 



Q = 2ttR Q_l 



Air'^RLjam-ySIiim-yS), (42) 



J = 2tiRJ^ + Jout =2TrR^m^SCbaka{f2 + fe) (43) 

and so we can eliminate fco in favor of J/luq and wq in 
favor of Q. Thus, the conserved quantities that enter into 
F always appear in the combinations Q and J / RQ. 
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At this point we have F as a function of the total 
charge Q = eZ, the angular momentum J, the vorton 
radius R, the string thickness S and the charged conden- 
sate at the center of the vortex \K~^\. We still need to 
minimize F in relation to R, S and while keeping 

Q and J fixed. Unfortunately, the expressions cannot be 
minimized analytically. We attack the question of the 
stabilization of the vorton numerically. We fix Z and J 
at some given values, and then numerically minimize E 
with respect to \K~^\, 6, and R simultaneously. 

For definiteness, we will show results for a gap A = 66 
MeV, a chemical potential /i = 450 MeV. We plot the 
resulting equilibrium radii in units of the string thickness, 
Ro/5 in Figure 3 as a function of Z and J. The photon 
mass normalized by the maximum possible photon mass 
m-y/ fe is shown in Figure 5, while and the total electrical 
current 1/ fe flowing in the vorton is shown Figure 4, all 
as functions of Z and J. 



40000 
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J 20 000 



10000 



Ro/S 




19 



10000 20000 30000 40000 



FIG. 3. The equilibrium radius Ro compared to the string 
thickness 5 at equilibrium as a function of charge Ze and 
angular momentum J. In all plots the region with Rq < S are 
shown in white. 



While the string thickness is a variational parameter 
and is varied for each (Z, J) pair, its change is very 
slight — for this choice of A and fi and quark masses, the 
value of S is varies from 25 to 29 fm. This small change 
is unimportant for understanding the presented figures — 
one may think of 5 as a fixed parameter. It is important 
to recognize, however, that 6 sets the scale for distance 
over which the fields vary appreciably. The numerical 
value we found is larger than the scale 1/A « 4 fm at 
which the effective theory in eq. (2) would break down. 
This justifies a posteriori the use of the effective theory 
in analyzing the physics of vortons. 
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FIG. 4. The current / flowing in a vorton of charge Ze and 
angular momentum J at its preferred radius. 
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FIG. 5. The photon mass for a vorton at its preferred radius 
with given charge Ze and angular momentum J. 



In Figure 3 we see that the equilibrium radius of a vor- 
ton increases with Z but is relatively independent of J. 
This is an artifact of choosing J as the vertical axis: the 
angular momentum always enters expressions divided by 
Z as noted above. However, these axes are useful for see- 
ing the other expected behavior. For example, in Figure 4 
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we see the current that flows through equiUbrium-sized 
vortons. Since at this scale, the radius is relatively con- 
stant with changing J, a higher angular momentum must 
directly imply a higher current. Such behavior is evident 
in Figure 4: for a fixed Z the current rises as J increases. 

From Figure 4 we see that larger Z lead to smaller cur- 
rents. On the other hand, larger Z correspond to larger 
i?o on account of Coulomb repulsion between opposite 
sides of the vorton (and also seen on Figure 3). Thus 
there is an inverse relation between current and radii: 
I ^ 1/R. This is the scaling one would naively expect if 
we assumed the current to be proportional to N / R as in 
the ungauged case. Qualitatively, this repulsive effect is 
still there in the more rigorous analysis presented here. 

To see the angular momentum barrier at work, we need 
to zoom out so that the suppression of J by Z is negligi- 
ble. To this effect, in Figure 6 we plot Ro/6 while fixing 
the variables Z and J/Z for the same choice of A and /i. 
One clearly sees the tendency for vortons to grow as J/Z 
is increased. 



Ro/d 




z 

FIG. 6. The equilibrium radius Ro compared to the string 
thickness 5 replotted with the vertical axis rescaled by Z. 

From this numerical example, typical among other rea- 
sonable values of the parameters, we learn some impor- 
tant lessons about vorton stability. The first and most 
important one is that in order to have a vorton that re- 
ally looks like a torus with R^ S charges of the order of 
Z many thousands are required. Second, even for these 
large values of Z the equilibrium radii are small enough 
that the energy from the membrane stretched across the 
vorton contributes negligibly to the energy budget. As 
mentioned before S seems large enough that the effective 
mesonic theory is believable. A more delicate issue is the 
charge neutrality of the medium as a whole. 



Given the requirement of overall charge neutrality for a 
physical system like a neutron star core, one expects there 
to be light charged leptons (either electrons or positrons) 
present in the bulk to cancel the electric charge of the 
vortons. Thus, electrons will orbit positively charged vor- 
tons. If the orbit size is much larger than the vorton size 
R, their orbits will have a typical size of the order of the 
Bohr radius aq divided by Z: R^ w ao/Z. For Z w 10000, 
the orbit size is sa 5 fm, smaller even than the thickness 5 
of the vorton. This means that electrons will likely shield 
appreciably the electric force helping to keep the vorton 
stable, a contribution not included in our calculation. 
On the other hand, only a fraction of the electron will be 
close enough to the vortons to shield the electric force. 
Just like in a regular atom, the exclusion principle pushes 
some of the electrons to much bigger orbits and they are 
very much less attracted to the vortons as the other elec- 
trons shield them from the effect of the vorton charge. 
One simple way to include some of the Pauli principle 
effects is to use the Thomas- Fermi approximation that, 
at least for a point nucleus, gives a typical electron orbit 
of the order of R^ « ao/Z^/^ « 2300 fm (for Z = 10000), 
which is much larger than the vorton radius i? w 120 fm. 
This suggests that only a small fraction of the electron 
orbiting a vorton will effectively shield the repulsion be- 
tween different parts of the vorton. The proximity of the 
scales discussed above, however, should serve as caution 
against hastily dropping the electron shielding effect. A 
better study of this problem would requires us to un- 
derstand the "atomic physics" of the vortons, where the 
toroidal shape of the "nucleus" and the electromagnetic 
fields sourced by the vorton have a significant effect on 
the electron orbitals. We will postpone this analysis to 
a future paper. If it turns out that the electron orbits 
are indeed much larger than the vortons, then "vortonic 
chemistry" will include all the molecules we would find 
on Earth if elements with Z « 10000 were stable. 

There is another class of potentially stable vortons 
which is not affected by the electron shielding. For small 
Z but very large angular momentum, say J « 10^, vor- 
tons can be stabilized at radii R ^ S. Since for smaller 
charges even the tightest electron orbits will be larger 
than the vorton and, in any case, electrons are less ef- 
fective at shielding magnetic forces, the electron cloud 
cannot destabilize the vortons. 



V. DISCUSSION 

We have shown that the CFL+if^ phase supports sta- 
ble charged current-carrying vortons, supported either by 
conservation of angular momentum, Coulomb repulsion, 
or a combination of both effects, and have estimated their 
size as a function of the relevant parameters as charge and 
angular momentum. Besides a more accurate calculation 
of the vorton structure, a number of improvements and 
future directions suggest themselves at this point. 

However, one crucial step is missing for connecting 
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all of these speculations with phenomenological conse- 
quences for neutron stars. Unfortunately, despite some 
effort, we were unable to understand and quantify the 
possible mechanisms of vorton formation during the early 
stages of neutron star evolution. For this reason we can 
not, at the moment, estimate their density. But all pos- 
sible phenomenological consequences of the existence of 
vortons on the physics of neutron stars depend crucially 
on the vorton density. In reference [17] a possible mech- 
anism for the generation of vortons was suggested, in- 
cluding the generation of a supercurrent. Still, no nu- 
merical estimate of the vorton density was made. Also, 
the crucial role of the charge in stabilizing the vorton 
was not understood at the time [18] and no mechanism 
for depositing charge along the vorton string is presently 
known. Still, if vortons are to exist in neutron star cores, 
it seems likely that they would have to be produced in 
non-equilibrium processes soon after the formation of the 
neutron star. If current and charge carrying vortons are 
indeed formed, the evidence we presented above strongly 
suggests they would be stable against strong and electro- 
magnetic processes. Estimating the number of vortons 
produced at the birth of the neutron star seems to us to 
be the most important problem to be addressed before 
further progress can be made in the topic. Another es- 
sential element missing in our analysis is an estimate of 
the vorton lifetime against weak processes, in particular 
the electron capture of inner electrons in a vortonic atom 



that can lower the charge and eventually destabilize the 
vorton itself. 

One may also imagine that a neutron star cored filled 
with vortons and electrons might possess crystalline or 
metallic structure due to the interactions mediated by 
the electron cloud or through magnetic dipole interac- 
tions. At the core of these vortonic atoms there is also 
a "nuclear structure" to be understood. We have ana- 
lyzed vortons in isolation and not considered either their 
interaction or their stability against fission into two vor- 
tons. Once a bulk with vortons is made, they should be 
locked in, supported by their conserved quantities and 
prevented from disappearing, except by weak effects and 
the possibility of fusing or fissioning. In the long term, 
then, one might expect all vortons to either decay or com- 
bine toward making the vortonic equivalent of ^^Fe. Most 
likely, the "chemical equilibrium" will never be reached, 
just like the chemical equilibrium among ordinary nuclei 
is never reached in the early universe. 
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